Anomalous reflection and excitation of surface waves in metamaterials 
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We consider reflection of electromagnetic waves from layered structures with various dielectric 
and magnetic properties, including metamaterials. Assuming periodic variations in the permittiv- 
ity, we flnd that the reflection is in general anomalous. In particular, we note that the specular 
reflection vanishes and that the incident energy is totally reflected in the backward direction, when 
the conditions for resonant excitation of leaking surface waves are fulflUed. 

PACS numbers: 42.25.-p 



Introduction Materials with new kinds of dielectric 
and magnetic properties, e.g. left-handed metamateri- 
als, have stimulated much recent work, see e.g. [ij. In 
the light of the anomalous electromagnetic characteristics 
[11, II 13, H, S B B it is necessary to reconsider 
many well known problems, accounting also for meta- 
materials. Due to their extraordinary properties, left- 
handed materials are believed to become highly impor- 
tant in various technological applications P, |3] ■ The first 
examples of metamaterials had their left-handed proper- 
ties in the microwave range, and then later for infrared 
waves . Recently materials showing anomalous refrac- 
tive behaviour have been constructed also in the optical 
regime Q. 

It has previously been shown in Ref. that certain 
layered dielectric structures give rise to peculiar reflection 
phenomena, provided surface waves 0,M E, 14, 111 can 
be excited. Such waves have previously been considered 
in connection with layered structures 0, [l^, [H, [Shopen 
plasma wave guides [13, [l3| , and fusion plasmas |l9| . 

In the present paper, we consider the reflective proper- 
ties of a dielectric and magnetic medium with arbitrary 
permittivity and permeability, which is covered by a di- 
electric layer. We show that for a homogeneous medium, 
we have the familiar specular reflection independently of 
the material parameters. However, with periodic varia- 
tions in the dielectric permittivity, as can be found in na- 
ture [2^, 21 1 and constructed in laboratories [H, d, [lH , we 
find that the incident energy can be completely converted 
into a backscattered energy flux, a phenomenon which is 
associated with the resonant excitation of leaking surface 
waves. The necessary conditions for obtaining backward 
reflection will be deduced, and possible applications of 
this peculiar mirror effect are consequently pointed out. 

Homogeneous media We first consider a semi-infinite 
{z > zq, region a) material with arbitrary relative per- 
mittivity Ea and permeability ^a- Specifically we will 
allow these constants to have arbitrary signs, such that 
for example metamaterials with negative signs of both e 
and ^ are included as special cases. This material is as- 
sumed to be covered by a dielectric layer with arbitrary 
permittivity and /ifc = 1 in the region < z < zq 
(region 6), whereas we have vacuum at z < 0. Fur- 



thermore, we let a p-polarized wave with a magnetic 
field {Boy /2)exp{ikxX + ik^z — iLut) + c.c, be incident 
on this structure. Here c.c. denotes complex conjugate, 
Bq is the wave amplitude, and kx.z > 0. We denote 
the angle of incidence by 9 (given by k^ ~ (cj/ c) sin 6 or 
kz — [ijJ / c) cos 9) . As a prerequisite for the more com- 
plicated calculations in the next section, we here study 
the fields induced in the materials for the special case 
where the waves are evanescent in both regions a and b. 
The magnetic field in each of the three regions z < 0, 
< z < Zq and z > zq is thus written in the form 
B = {y /2)B(x, z)exp(— iwi) -I- c.c. In the vacuum region 
the y-component satisfies the equation 



d^B + dlB 
which has solutions 



-B = 0, 



(1) 



B = Bne 



Bo {Rs 



(2) 

where Rs and Rf, are constants which represent the spec- 
ularly and backward reflected waves, respectively. Nor- 
mally, one would expect i?^ to be zero. However, we note 
that metamaterials may behave differently compared to 
normal materials. Furthermore, inclusion of the i?b-term 
will be useful for the calculations made in the next sec- 
tion. From energy conservation, evanescent waves in re- 
gions a and h imply Rb + Rs = 1. Next, in region b we 
solve the equation 



dlB + djB+'^B^Q, 



and write the solution in the form 



(3) 



2 \ Kb J 



+ [Rse' 



(4) 



where Kb = {lo / c){s\t? 9 — et)^/^, and where r = {\ + 
ikzEb/ Kb) / (\ — ikzEb/Kb) represents the reflectivity at 
z — 0. We note that the standard boundary conditions 
expressed in terms of B, namely the continuity of B and 



2 



{l/e)dzB are satisfied at z = for arbitrary values of Rg 
and Rh. In region a we use the equation 



d^B + dlB - 
and write the solution as 



B^O, 



(5) 



(6) 



with Ka = (a;/c)(sin^0 — Ea/ia)^/^. Since both /i and e 
vary across the boundary, the boundary conditions are 
now the continuity of B/n and (l/e)dz{B/ n). Matching 
the solutions in regions a and b, we obtain 



2 V 



B4 



+ (e'^f^o + re-'^^^") i?J = (7) 

Ma 



Kb 



-^''°)Rb = —, (8) 

Ma 



as well as 



Kb 



+ (e" 



bZO „c~''''^Ol 



(9) 



Kb 



-Kb.o)^^^„^^i£i^lO) 

^aMa 



Here we note that the equations involving B^ couple only 
to the backward reflection coefficient Rb, and that the 
incident wave does not act as a source for this part. Thus, 
we can interpret this contribution as a leaking surface 
23 



wave 



22 



that is localized around the boundary z = zq 
and decays exponentially away from it, but couples to 
a propagating ("leaking") part in the boundary region. 
However, no matter what kind of material properties we 
have in region a, this mode cannot be excited by the 
incident wave, as the boundary conditions do not allow 
it to couple to the incident wave. The field profile of the 
surface mode is shown in Fig. 1. Furthermore, as will be 
demonstrated in the next section, the properties of the 
surface wave are crucial for the reflection properties when 
the media in the different regions are not homogeneous. 

Due to the decoupling of the surface waves in a ho- 
mogeneous medium, however, we here obtain Rg = 1 
from the boundary conditions, and hence we always have 
a purely specular reflection with this geometry. Let us 
stress that for a metamaterial, the surface wave localized 
around z = zo carries energy in the opposite direction 
along the x-axis, as compared to the incident wave. How- 
ever, as shown above, this does not affect the direction 
of the reflected energy in the vacuum region. 

Periodic media We next extend our previous analysis 
to consider a semi-infinite (z > zq) medium with a pe- 
riodically varying permittivity Ea = £aO + £ai cos{2kxx), 




z=0 



FIG. 1: A schematic representation of the z-dependence of 
the leaking surface wave in the various regions. 



where Eao and £ai are constants. We note that such a 
medium can be constructed in laboratories [l|, [^ .[lll , but 
also that such structures are common in nature [20| . The 
medium in region b is the same as in the previous section. 
For simplicity, the permeability jia in region a is kept con- 
stant. Generally, the interaction between the incident 
wave and the periodic dielectric structure in medium a 
will result in the generation of new waves corresponding 
to the sum and difference of the wave vectors. Here, we 
choose the a;-component of the incident wave vector such 
that it is equal to half the wave number of the periodic 
modulation in region a. This choice of incident angle 
9 means that the backscatter generation mechanism is 
resonant and thus much more important than the cor- 
responding interaction mechanism for other values of 9. 
Since the medium and thereby the wave equations are 
the same as in the previous analysis in the vacuum region 
and region b, we can adopt the corresponding solutions 
(jll) and ([6]). Furthermore, for region a we can closely 
follow the analysis in Ref. 12]. Including the spatial 
dependence of the permittivity, the wave equation ^ is 
then generalized to 



-^dlB + dJ-^dxB] + '^B = Q. (11) 

Similarly, the ansatz ^ is generalized to 

B = B+{z)e''''''' + B_(z)e-*'===^. (12) 

Assuming that Sai is much smaller than £^01 we then 
obtain the two coupled equations 



dtB± ^ Ki~ 



, ,2 



B± 



2£aOMa 

with the boundary conditions 



[S'^T - {kl + 2^1) St] , (13) 



{B- 



+ )z=za 



Bo [r + Rs + il + rRs 



-2KbZol 



(14) 
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, Sal 



Rh I 



BoRb (l - 
'o — (eaoMa [r 

Eh 



(15) 



Rs 



(1 + ri?,)e^2K6^o] 
(16) 



and 



{dzB^)z=za = Ba — leaOf^aRb (l - 

+^ [r + i?.-(l + ri?,)e-2«.^o)] 1^ (17) 

where Ka = (w/c)(sm^ 6'-eaoMa)^^^ and Bq = (-So/2)(l- 
jfc2£6/Kb)exp(Kf,0o) The solution of (|13|) can be found by 
expanding Bj- in powers of the small parameter Sai ■ Af- 
ter lengthy but straightforward calculations, we obtain 
the solution 



B± w C±e- 



(18) 

where C±, as well as Rb and i?^, are complex constants 
that can be found from (fT4|) -p7 l) . As Sai/Sao ^ 1, only 
weakly damped leaking surface waves are of interest here. 
Thus, from now on we focus our attention on the regime 
exp(— KdZo) <C 1. We then obtain the simple approxi- 
mate solution 



Rs 



£12 



1 



2^ -2k(,zo 



and 



where 



Rb 



D 



Eal 



co^ (l-r2) 



2k2c2 



D2 



l^aEb 



(19) 



(20) 



(21) 



In the absence of dielectric modulations, i.e. if Sai = 
0, we recover from (jl9l) the well known formula = 
-r|D|2/D2, i.e. \Rf\ = 1. Obviously, = for this 
case. However, from Eqs. (fT9| and ([20|) . we also find the 
interesting result that i?s ~ and Rb ~ ir if the real 
part of the surface wave dispersion function ([2T|l is equal 
to zero, i.e. 



1 

and if 



Kb£aO^J■a 



2(l^iMAi)^-2K,.o 



KaSb {l + klel/nD 

8^2 c2 



0, (22) 



Kb{l + klel/Kl) 



(23) 



When i?s ~ it also follows that w 1. Equation 
means that the incident wave resonantly excites a leaking 



surface wave, in which most of the energy is concentrated 
near the interface at z = zq. Equation (|23p then defines 
the condition for the backward energy flux to be equal 
to the incident flux. We note that the relation only 
can be fulfilled if eaoHa/sb < 0. For example, for a meta- 
material in region a {Sao < 0, < 0) we must have neg- 
ative dielectric permittivity also in region 5, whereas for 
a nonmagnetic material the dielectric permittivity must 
change sign between the regions a and b. In general, to be 
able to study the present phenomenon, it is necessary to 
consider configurations with suitable properties so that 
weakly damped leaking waves can exist [lj| . Specifically 
in the absence of the boundary region (i.e. if zq = 0, or 
if £(, ~^ 1); there are no such waves. In our model, it is 
possible to investigate leaking waves, however [l3|. Simi- 
lar studies can naturally be performed for other bounded 
media, but the mathematics would then be much more 
involved. 

Summary In the present paper, we have considered 
the electromagnetic reflection properties of an arbitrary 
dielectric and magnetic material covered by a dielectric 
layer. For homogeneous media, it was shown that leak- 
ing surface waves cannot be excited by an incident wave, 
and we have thus the familiar specular reflection. For 
small periodic variations of the dielectric permittivity 
^2^1, however, the situation is radically different, and the 
coupling to leaking surface waves can then strongly in- 
fluence the reflective properties. In particular, we have 
found that specular reflection is absent {Rs = 0) if the 
frequency of the incident wave satisfies the dispersion re- 
lation ((22l) . and if the width zq of the dielectric layer is 
related to the amplitude Sai of the permittivity modula- 
tion by means of the relation . All the incident wave 
energy is then reflected in the backward direction. These 
properties make it possible to pick out a specific angle 
and a specific frequency for a wave which is backward 
reflected, while still keeping most of the signal in the or- 
dinary reflected signal. Obviously, such a peculiar mirror 
effect is of experimental interest. 



[1] 

[2] 
[3] 

[4] 

[5] 

[6] 

[7] 



V. M. Agranovich and Yu N. Gartstein, Physics-Uspekhi, 
49, 1029 (2006). 

V. M. Shalaev, Nature Photonics, 1, 41 (2007). 

V. P. Makarov and A. A. Rukhadze, JETP, 103, 354 

(2006). 

A. Chowdhury and J. A. Tataronis, Phys. Rev E, 75, 
016603 (2007). 

M. Marklund, P. K. Shukla, L. Stenflo, and G. Brodin, 
Phys. Lett. A, 341, 231 (2005). 

M. Marklund, P. K. Shukla, and L. Stenflo, Phys. Rev. 

E, 73, 037601 (2006). 

R. Ruppin, Phys. Lett. A 277, 61 (2000). 

I. V. Shadrivov, A. A. Sukhorukov, Yu. S. Kivshar, A. 

A. Zharov, A. D. Boardman, and P. Egan, Phys. Rev. E 

69, 016617 (2004). 



4 



[9] I. V. Shadrivov, A. A. Zharov, and Yu. S. Kivshar, Appl. 

Phys. Lett. 83, 2713 (2003). 
[10] I. V. Shadrivov, R. W. Ziolkowski, A. A. Zharov, and Yu. 

S. Kivshar, Opt. Express 13, 481 (2005). 
[11] K. Sanderson, Nature 445, 346 (2007). 
[12] Yu. M. AHev, A. A. Frolov, G. Brodin and L. Stenflo, 

Phys. Rev. E, 47, 4623 (1993). 
[13] O. M. Gradov and L. Stenflo, Phys. Scripta, 25, 631 

(1982). 

[14] O. M. Gradov and L. Stenflo, Phys. Rep. 94, 111 (1983). 

[15] L. Stenflo, Phys. Scripta T63, 59 (1996). 

[16] I. Zhelyazkov, The Physics of Ionized Gases, edited by J. 

Puric and D. BeUc (World Scientific, 1987), p. 459. 
[17] D. Grozev, A. Shivaxova and S. Tanev, J. Plasma Phys. 



45, 297 (1991). 

[18] Yu. M. Aliev, H. Schliitcr and A. Shivarova, Guided wave 
produced plasmas, (Springer, Borhn 2000). 

[19] M. Y. Yu, Surface Waves in Plasmas and Solids, edited 
by S. Vukovic (World Scientific, 1986), p. 419. 

[20] P. Vukusic, Phys. Today, 59, 82 (2006). 

[21] L. Stenflo and I. P. Yakimenko, Phys. Scripta 18, 151 
(1978). 

[22] A. V. Kats, S. Savel'ev, V. A. Yampol'skii, and F. Nori, 

Phys. Rev. Lett. 98, 073901 (2007). 
[23] E. Fourkal, I. Velchev, C.-M. Ma, and A. Smolyakov, 

Phys. Lett. A 361, 277 (2007). 



